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ABSTRACT 


The  effects  of  Gausslan-shaped  IF  and  video  filters  on  the 
performance  of  a  square-law  detector  employing  post-detection 
integration  are  analyzed.  The  number  of  additional  noise-only 
samples  that  are  integrated  due  to  a  finite  video  bandwidth  is 
determined.  Emerson's  method  is  used  to  obtain  an  expression 
for  the  cumulants  of  the  output  probability  density  function. 
These  cumulants  are  used  in  Edgeworth's  asymptotic  series 
expansion  of  the  density  functions.  By  integrating  these 
density  functions,  the  receiver  operating  characteristics  are 
determined  for  various  ratios  of  IF  bandwidth  to  video 
bandwidth.  A  collapsing  ratio  which  takes  into  account  the 
Gaussian  shape  of  the  filters  is  formulated  and  compared 
against  Barton's  approximate  formula  for  the  collapsing  ratio. 
For  typical  video  bandwidths,  Barton's  approximate  formula  is 
found  to  overestimate  the  collapsing  loss  by  an  amount  less 
than  0.5  dB. 
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I. 


INTRODUCTION 


This  report  is  concerned  with  the  effects  of  insufficient 
video  bandwidth  on  the  detection  performance  of  a  square  law 
detector  with  Gaussian-shaped  intermediate  frequency  (IF)  and 
video  filters  in  pulsed  radar  applications. 

The  square- law  detector  is  assumed  to  have  an  arbitrary 
wide  low-pass  zonal  filter  [Ref.  1]  so  that  its  instantaneous 
output  is  just  the  square  of  the  envelope  of  the  input  signal. 
This  assumption  is  in  line  with  Marcum's  definition  of  the 
detector  as  "any  device  whose  instantaneous  output  is  a 
function  of  the  envelope  of  the  input  wave  only”  [Ref.  2]. 

In  most  practical  radars  with  superheterodyne  receivers, 
the  second  detector  is  usually  followed  by  one  or  more  stages 
of  video  amplification  prior  to  post-detection  integration  and 
threshold  detection.  The  bandwidth  of  the  video  amplifier 
must  be  sufficiently  wide  to  avoid  distortion  of  the  input 
video  signal  from  the  detector.  It  has  been  pointed  out  that 
the  video  bandwidth  By  must  be  greater  than  one-half  the  IF 
bandwidth  to  be  "sufficiently  wide"  [Ref.  3].  In 
practice,  the  video  bandwidth  may  be  less  than  half  the  IF 
bandwidth . 

An  insufficient  video  bandwidth  will  have  two  effects  on 
the  detection  statistics  of  the  radar  receiver.  Firstly,  it 
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will  affect  the  threshold  setting  for  a  specified  probability 


of  false  alarm,  Secondly,  it  will  adversely  affect  the 
input  signal-to-noise  ratio  required  to  attain  a  desired 
probability  of  detection,  P^,  for  a  specified  probability  of 
false  alarm.  An  equivalent  statement  to  the  second  effect  is 
that  there  will  be  a  degradation  in  probability  of  detection 
vis-a-vis  the  probability  of  detection  attainable  with  an 
ideal  video  bandwidth  for  the  same  signal-to-noise  ratio  and 
probability  of  false  alarm. 

The  focus  of  this  report  will  be  on  the  degradation  in 
performance  associated  with  the  second  effect.  Mathematically, 
Marcum  [Ref.  2]  has  treated  this  degradation  in  performance  by 
modelling  it  as  a  consequence  of  integrating  M  additional 
noise-only  samples  along  with  the  desired  N  signal-plus-noise 
samples^.  The  degradation  in  performance  is  quantified  by 
introducing  a  "collapsing  loss".  It,  which  has  been  defined  in 
dB  as  the  incremental  increase  in  signal-to-noise  ratio  with 
respect  to  the  signal-to-noise  ratio  required  in  the  ideal 
case  of  integrating  only  signal-plus-noise  samples  for  the 


’^Marcum's  paper  first  introduced  the  general  concept  of  a 
collapsing  loss  resulting  from  the  integration  of  a  greater 
number  of  noise-only  samples  than  signal-plus-noise  samples. 
Insufficient  video  bandwidth  is  treated,  among  others,  as  a 
specific  cause  of  collapsing  loss. 
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same  probability  of  detection.  Mathematically,  this  is 

expressed  as:  [Ref.  4] 

Lp  dB  -  dB  -  Xi^dB  ( 1 ) 


where  X„+h  is  the  signal-to-noise  ratio  required  when  N  signal- 
plus-noise  samples  plus  M  extra  noise  samples  are  integrated, 
and  X„  is  the  signal-to-noise  ratio  required  when  only  N 
signal-plus-noise  samples  are  integrated. 

The  signal-plus-noise  samples  and  noise-only  samples  are 
assumed  to  be  independent  and  respectively  identically 
distributed. 

Mathematically,  Marcum  showed  that  for  a  non-fluctuating 
target,  the  effect  on  the  detection  statistics  of  integrating 
M  additional  noise  samples  can  be  taken  into  account  by  multi¬ 
plying  the  joint  characteristic  function  of  the  N  signal-plus- 
noise  samples  by  the  joint  characteristic  function  of  the  M 
noise-only  samples.  Defining  the  collapsing  ratio, 


P  - 


M  *  N 
N 


(2) 


Marcum  showed  that  for  the  square  law  detector 


(5)  - 


exp[  -Np  (  — )  ( — — rg- )  3 


(3) 
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where  Ch*k(  € )  is  the  characteristic  function  for  the  joint 
density  function  of  N  signal-plus-noise  samples  plus  M 
additional  noise  samples. 

When  p  is  equal  to  one  (i.e.,  M  is  equal  to  zero),  Eqn 
(3)  reduces  to  the  equation  of  the  characteristic  function 
with  no  collapsing  loss,  i.e.. 


exp[  -NX(-^)  ] 
_ 1+J^ 


/ 1  .  c  \  W 


Comparing  Eqn  (3)  with  Eqn  (4),  it  is  seen  that  the 
statistical  results  for  p  =  1  can  be  easily  extended  for  any 
p.  In  particular,  receiver  operating  characteristics  (ROC) 
such  as  Meyer  plots^  for  p  =  1  are  available  in  the 
literature  and  can  be  used  directly  to  determine  the  signal- 


to-noise  required  by  using  pN  instead  of  N  as  the  number  of 


pulses  integrated  and  with 


—  as  a  parameter  of  the  ROC 


instead  of  X.  [Ref.  4] 


^These  are  plots  of  as  a  function  of  the  number  of 
pulses  noncoherently  integrated  for  a  given  Fehlner's  false 
alarm  number  and  with  the  input  signal-to-noise  ratio  for  a 
single  pulse  as  a  running  parameter. 
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Hence,  the  quantitative  effects  on  the  detection 
statistics  of  insufficient  video  bandwidth  can  be  determined 
once  the  collapsing  ratio,  p,  is  known.  The  problem  is  in 
relating  M  (and  hence,  p)  to  the  video  bandwidth.  What  is 
the  effective  number  of  additional  noise-only  samples  that  are 
integrated  along  with  the  intended  number  of  signal-plus-noise 
samples  given  a  finite  video  bandwidth? 

A  common  practice  is  to  use  the  following  semi-empirical 
formula  given  by  Barton  [Ref.  5]  to  estimate  the  collapsing 
ratio: 


Pb 


Bjf  ^  2B, 
2B„ 


(5) 


Where  the  IF  noise  equivalent  bandwidth  Bi,  is  assumed  to  be 
matched  to  the  pulse  width  such  that  Bi,t  =  l.  The  video 
filter  is  modeled  as  an  ideal  low-pass  filter  with  a  flat 
frequency  response  from  -B,  to  +Bv  and  zero  elsewhere.  The 
frequency  response  characteristics  of  the  actual  filter  is  not 
taken  into  account  apart  from  computing  its  noise-equivalent 
bandwidth,  Bv. 

The  purpose  of  this  paper  is  to  investigate  the  effects  of 
Gaussian- shaped  IF  and  video  filters  on  the  detection 
performance  of  the  square-law  detector.  The  primary 
consideration  in  the  choice  of  Gauss ian-shaped  filters  for 
this  investigation  is  the  mathematical  convenience  of  the 
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Gaussian  shape  which  yields  analytical  results  in  closed  form. 

Although  a  Gaussian  filter  is  physically  not  realizable, 
it  has  been  found  to  be  a  good  approximation  in  characterizing 
the  frequency  response  of  multistage  RC  amplifiers  [Ref.  6]. 
Hence,  the  results  obtained  will  have  direct  application  for 
square-law  detectors  with  multistage  pre-detection  and  post¬ 
detection  RC  filters,  besides  providing  a  relative  measure  of 
the  accuracy  of  Barton's  formula  for  the  collapsing  ratio. 

In  the  next  section,  the  effects  of  video  bandwidth  on  the 
detector's  output  probability  density  functions  and  operating 
characteristics  will  be  analyzed.  Specifically,  a  formula 
will  be  developed  to  determine  the  additional  noise-only 
samples  integrated  due  to  a  finite  video  bandwidth.  This 
leads  to  the  formulation  of  the  collapsing  ratio  which  will  be 
discussed  and  compared  against  Barton's  formula  in  Section 
III.  Section  IV  summarizes  the  methods  and  findings  of  this 
paper . 
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II.  DEGRADATION  IN  PERFORMANCE 


A.  DESCRIPTION  OF  THE  SYSTEM 

The  system  under  consideration  is  illustrated  in  Figure  1. 
It  consists  of  a  Gaussian  IF  filter,  a  square-law  detector,  a 
video  amplifier  with  a  Gaussian  frequency  response,  and  a 
linear  post-detection  integrator.  For  convenience,  it  shall 
be  referred  to  as  the  Gaussian  square-law  system. 


Figure  1.  Block  diagram  of  the  Gaussian  square-law  system. 


The  IF  filter  has  a  Gaussian  passband  characteristic  with 
a  center  frequency  of  and  a  root-mean-square  bandwidth  of 
P  .  The  video  amplifier  also  has  a  Gaussian  passband  but  with 
a  root-mean-square  bandwidth  of  v  centered  at  zero  frequency. 
Since  the  gain  of  the  video  amplifier  will  have  no  impact  on 
detection  performance,  it  is  assumed  to  have  unity  gain  for 
mathematical  convenience.  Correspondingly,  the  impulse 
response  functions  of  the  IF  filter  and  video  amplifier  are 
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respectively  defined  by: 


— 

h^fit)  -  2(2np2)2  exp[- (2iip)2_]  cosw^t 

and 

— 

h^it)  -  (27iv2)  2  exp  [  -  (2itv)2  _  j  (7) 

The  signal  to  be  detected  is  a  pulse-modulated  sinusoidal 
waveform  corresponding  to  the  echo  from  a  target  at  a 
particular  range  from  the  radar. 

The  input  noise  to  the  system  is  assumed  to  be  additive 
white  Gaussian  noise  (AWGN)  with  a  two-sided  power  spectral 
density  equal  to  x\/2>  The  noise-equivalent  bandwidth  of  the 
IF  filter,  Bit,  is  matched  to  the  pulse  width  so  that  noise  at 
the  input  to  the  detector  can  be  assumed  to  be  completely 
uncorrelated  from  pulse  to  pulse  [Ref.  2]. 

B.  ADDITIONAL  NOISE-ONLY  SAMPLES  INTEGRATED 

For  the  hypothetical  case  of  infinite  video  bandwidth,  the 
video  range  bins  will  have  the  same  resolution  as  the  IF  range 
bins,  which  is  of  the  order  of  one  pulse  width.  However,  for 
most  practical  video  amplifiers  with  finite  bandwidths,  the 
resolution  of  the  video  range  bins  will  be  poorer  than  those 
of  the  IF  range  bins.  In  such  cases,  each  sampled  output  of 
the  video  amplifier  will  be  shown  in  the  subsequent  paragraphs 


8 


to  be  equivalent  to  the  result  of  averaging  independent 
samples  of  the  detector  output,  where  is  a  number  greater 
than  one.  Since  the  detector  follows  a  square  law,  each  of 
the  Nv  independent  samples  is  equal  to  the  square  of  the 
voltage  envelope  in  the  corresponding  IF  range  bins. 

Suppose  a  target  is  present  at  a  particular  IF  range  bin. 
One  of  the  Nv  samples  will  consist  of  signal-plus-noise.  The 
remaining  Nv-1  samples  will  consist  of  noise  only.  Hence,  Nv“l 
additional  noise-only  samples  will  be  integrated  for  every 
signal-plus-noise  sample  integrated.  The  added  noise  results 
in  a  degradation  in  performance  that  is  quantified  by  a 
collapsing  loss. 

The  number  of  noise-only  samples  integrated  along  with  the 
signal-plus-noise  sample  can  be  determined  by  recognizing  that 
the  video  amplifier  with  a  finite  bandwidth  behaves  as  an 
analog  filter  integrator  and  paralleling  it  with  a  digital 
integrator  that  averages  Nv  independent  samples  of  the 
detector  output.  This  is  illustrated  in  Figure  2. 

The  digital  integrator  as  shown  in  Figure  2(a)  has  the 
well-known  effect  of  decreasing  the  variance  of  the  input 
samples,  q^,  by  a  factor  of  Nv,  i.e.. 


var  {y) 


Var  {<7j} 
K 


(8) 


where  y  is  the  output  of  the  digital  integrator. 
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(a)  Digital  integrator:  \/ar(y)  •  Var(q)/N 


(b)  Video  filter:  Equivalent  -  Var{q(t))/Var{y(t)) 


Figure  2 .  Equivalent  number  of  independent  samples  averaged 
by  the  video  filter. 


Hence,  the  equivalent  number  of  samples,  Nv,  averaged  by 
the  video  filter  can  be  determined  by  evaluating: 


Var  {g{t) 

Var  ly  ( t) }  ( ^ 


where 

q(t)  =  output  of  the  square-law  detector 
y(t)  =  output  of  the  video  filter  (amplifier). 

The  cumulants  of  the  probability  density  function  (pdf)  fy 
at  the  output  of  the  Gaussian  video  filter  can  be  determined 
by  the  method  of  Emerson  (Ref.  7].  When  the  input  to  the 
square  law  detector  is  a  CW  sinusoidal  signal  contaminated  by 
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AWGN,  the  cumulants  are  given  by: 


K. 


•y-r 


_ (2a^)  (r-1)  » _ 

j, 

[  (1  +  27^)  2  +  1]  ^  -  [(1  +  2y2)  2  _  1]  r 


(10) 


r 

L 


1 


+  rX' 


[  (1  +  2y^)  ^  1]  ^  - 

[  (1  +  2y2)  2  +  1]  r  + 


[(1-^2y^)  ^ 
[(142y")  " 


il! 

1]  ^ 


_1 

(1+2y^)  ^ 


where 

=  rth  cumulant  of  the  pdf,  fy 
=  average  noise  power  at  the  input  to  the  detector 
X'  =  signal-to-noise  ratio  at  the  input  to  the  detector 


The  mathematical  derivation  of  Eqn  (10)  is  shown  in 
Appendix  A.  For  the  purpose  of  this  Section,  it  is  sufficient 
to  point  out  that  the  derivation  takes  into  account  the  shape 
of  the  Gaussian  filter  and  that  the  cumulants  are  in  general 
related  to  the  characteristic  function  of  the  associated 
probability  density  function  by  the  following  definition 
[Ref.  8]: 


exp[  ^ 


r-1 


K,(jV  ^  1 

r!  J 


-  C(5) 


(11) 
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The  curoulants  are  also  related  to  the  central  moments  of 


the  density  function.^  In  particular,  the  variance  or  second 
central  moment  of  the  density  function  is  equal  to  the  second 
cumulant  of  the  density  function. 

Hence,  for  the  case  of  a  CW  sinusoidal  signal. 


Variyit))  -  Ky^ 


y/l  +  2y^ 


r 

I 


1  +  2X' 


1 

1  +  i 


(12) 


To  use  Eqn  (12)  for  evaluating  in  Eqn  (9),  X'  must  be 
the  equivalent  CW  signal-to-noise  ratio  that  will  produce  the 
same  variance  resulting  from  the  averaging  of  one  signal-plus- 
noise  sample  and  N^-l  noise-only  samples. 

Referring  to  Figure  2(a),  let  X  be  the  input  signal-to- 
noise  ratio  of  the  jth  range  bin  containing  the  target.  Then, 


y  - 


i  *  i-  £ 


(13) 


and 


Var  {y) 


Var  {qr^.}  +  Var  (g,)  , 

nI  nI 


^  *  j  (14) 


’A  full  treatment  of  the  relationship  between  the 
cumulants  and  the  moments  of  the  probability  density  function 
can  be  found  in  Ref.  8. 
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where  are  independent  and  identically  distributed  noise- 
only  samples  for  i  =  j. 

Now,  is  the  output  of  the  square-law  detector  prior  to 
any  video  filtering.  For  analysis,  this  is  theoretically 
equal  to  the  output  y  when  the  video  filter  has  an  infinite 
bandwidth .  Hence , 


Var  {qj} 


^y-0,x'-x 


-  a^(l  +  2X) 


(15) 


and 


Var  {qj} 


I  Y-o,  x'-o 


(16) 


a 


2 


Substituting  the  above  results  into  Eqn  (14)  and  simplify 
ing  yields  the  following  expression  for  the  variance  of  y; 


Var  {y) 


X  + 


(17) 


Let  X  be  the  equivalent  signal-noise-ratio  in  each  of  the 
Wv  range  bins  that  will  produce  the  same  variance  resulting 
from  averaging  the  single  signal-plus-noise  sample  with  Wv 
noise-only  samples. 
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Then,  using  Eqn  (12)  with  y  “  0  and  X' •  X,  it  follows  that 


vax  [q^]  -  {1-^2X) 


and 


Ar„ 


Var  {y}  - 


li  E  -Ji  1 

■^''v  i-1 


-  —  (1  +  2^:) 


(18) 


(19) 


By  equating  Eqn  (17)  and  Eqn  (19),  it  can  be  easily  shown 

that 


X  - 

N„ 


(20) 


Hence,  for  an  arbitrary  y,  it  follows  from  Eqn  (12)  and 
Eqn  (20)  that 


Variyit)}  -  -^^,2  | 


[ 


(21) 


X 


_  ,  1  +  2-^., 

y/l  +  2y^ 


1  +  2y' 
1  + 
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The  equivalent  number  of  independent  samples  averaged 
by  the  video  filter  can  now  be  determined  by  using  Eqn  (9), 
i  .e.  / 


Var  (g(  t) 

Var  iy(t) }  (22) 


with 


Var  {y  ( t) )  - 


v/TTIr 


N  *  o  ^ 
1  +  2 


1  f2Y^ 


1  +  Y‘ 


(23) 


and 


V-ar{g(t)}  - 


N„ 


(24) 


Substituting  Eqn  (23)  and  Eqn  (24)  into  Eqn  (22)  results 
in  the  following  quadratic  expression  in 


[ 


N‘  +  i  2X, 


Y-i  -  ^/T72^]n^  -  2X)/1  +  2y^  -  0  (25) 


l+Y*^ 


Hence,  the  equivalent  number  of  independent  samples 
averaged  by  the  Gaussian-shaped  video  filter  is  given  by; 


K  - 


-b  +  yJb^~AC 


(26) 
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where 


b 


and 

c 


2X 


M 


l+2y^ 
1  + 


-2X  y/l  +  2y^ 


(27) 


(28) 


Solving  Eqn  (26)  for  y-0  shows  that  is  equal  to  one 
when  the  video  bandwidth  is  infinitely  large.  This  result  is 
consistent  with  the  earlier  observation  that  the  video  range 
bins  will  have  the  same  resolution  as  the  IF  range  bins  when 
the  video  bandwidth  is  infinite. 


C.  RECEIVER  OPERATING  CHARACTERISTICS 

The  collapsing  loss  caused  by  insufficient  video  bandwidth 
can  be  observed  by  examining  the  effect  of  integrating  the 
additional  noise-only  samples  on  the  output  probability 
density  function  and  operating  characteristics  of  the  Gaussian 
square  law  system. 

The  number  of  signal-plus-noise  samples  N  integrated  by 
the  post-detection  integrator  corresponds  to  the  number  of 
pulses  returned  from  a  target  as  the  radar  antenna  scans 
through  its  beamwidth  [Ref.  3].  For  most  practical  radars, 
this  number  is  sufficiently  large  so  that  the  probability 
density  function  of  the  voltage  at  the  output  of  the  linear 
post-detection  integrator  shown  in  Figure  1  tends  towards  a 
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normal  distribution  which  can  be  asymptotically  expanded  using 
Edgeworth's  groupings  of  the  Gram-Charlier  series  expansion 
[Ref.  9,  Ref.  10]: 


fz{z)  “  [4>(z') 


(29) 


+  ±1  (2/)  +  A  4,(6)  {z>) 

24  ^  72  ^ 


^  (z')  -  (z')  -  (z') 


120 


144 


1296 


where  the  terms  on  each  line  within  the  brackets  are  of  the 
same  order  of  magnitude  [Ref.  11];  and 


z  - 


z'  - 


integrator  output  voltage  normalized  by  the  rms 
noise  voltage  a  at  the  input  to  the  detector. 


oz  -  NK, 


■y.i 


^  (z)  - 


1  r  z‘ 

exp  p  — 


^/2‘n 


(30) 


(31) 
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^<m)  (2)  .  (-1)"  <|)(0)  i4(z) 

H^(z)  =  mth  order  Hermite  polynomial  in  z 


1 

NO.  Sr 


V,  Z*2 

)l.o.5r 


(32) 


The  symbol  denotes  the  rth  cumulant  of  the 

probability  density  function  of  the  video  filter  output  y 
which  takes  into  account  the  additional  noise-only  samples 
averaged  along  with  each  signal-plus-noise  sample. 

Studies  by  Cramer  [Ref.  11]  have  validated  the  asymptotic 
behaviour  of  Edgeworth's  series  and  determined  that  in 
general,  the  remainder  term  of  the  expansion  is  of  the  same 
order  as  the  first  term  neglected.  Hence,  the  output 
probability  density  function  fz(z)  can  be  expanded  to  the 
desired  order  of  accuracy  once  the  cumulants  JCy ,,  are  known. 

The  cumulants  Ky,^  can  be  determined  from  Eqn  (10)  with  X' 

X 

equal  to  — ,  i.e.. 
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Once  the  probability  density  function  f:(z)  is  determined, 
the  cumulative  distribution  function  Fi(Zy)  can  be  found  by 
integrating  fzCZiJ  in  Eqn  (30)  from  -<»  to  z^.  Hence, 


Fz(Zy)  -  1  P(z') 


(34) 


+  A  (j,'3)  (2/)  +  £l  A'B)  (2O 

24  ^  72  ^ 


■y  C  C 

4>(4)(2/)  _  4>(6)(z0  -  (2/) 

^  '  “i  A  A  ^  '  ' 


120 


144 


1296 


Where 


_  q  Zy  -  iVFy^  1 


yfNK^: 


(35) 


P{zO  -  p  — ^  exp  I  - 
v/ZiT  ^ 


f  z^  1 


j  dz  -  0 . 5  erf 


r  z'  ’ 


L 


v/2 


(36) 


The  probability  of  false  alarm  P,.  and  the  probability  of 
detection  P^  for  a  normalized  threshold  voltage  of  z^  can  then 
be  determined  by  taking  the  complement  of  Eqn  (34),  i.e.. 


•Pfa  “  ^  ^z,  0  ^ 


(37) 


Pd  -  1  - 


(38) 
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where  denotes  the  cumulative  distribution  for  the 
hypothesis  that  the  target  is  absent  while  F2,i(z^)  denotes  the 
cumulative  distribution  for  the  hypothesis  that  the  target  is 
present . 

The  receiver  operating  characteristics  is  then  obtained  by 
plotting  the  probability  of  detection  versus  the  probability 
of  false  alarm  with  X,  y,  and  N  as  parameters. 

It  is  obvious  from  the  above  expressions  that  the  process 
involved  in  the  computation  of  the  output  probability  density 
function  f(z)  and  the  receiver  operating  characteristics  is 
numerically  tedious.  However,  the  steps  involved  are  straight 
forward  and  can  be  easily  programmed  for  execution  in  a 
digital  computer.  As  part  of  this  thesis,  routines  in  PC- 
MATLAB  have  been  written  to  compute  the  cumulants,  the 
probability  density  functions,  and  the  associated  receiver 
operating  characteristics.  These  routines  are  listed  in 
Appendix  B. 

Using  these  routines,  the  effect  of  insufficient  video 
bandwdith  on  the  the  detection  statistics  of  the  Gaussian 
square-law  system  can  be  investigated  for  various  values  of  y , 
input  signal-to-noise  ratio,  and  W.  This  is  illustrated  in 
Figure  3  to  Figure  5. 

Figure  3  compares  the  output  probability  density  functions 
obtained  for  a  finite  video  bandwidth  equal  to  a  quarter  of 
the  IF  bandwidth  (y~2)  against  the  corresponding  density 
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N-9 


Figure  3.  Gaussian  square  law  system  output  probability 
density  function. 

functions  when  the  video  bandwidth  is  infinite  (Y“0).  The 
output  probability  density  functions  corresponding  to  the 
absence  and  presence  of  a  target  are  respectively  denoted  by 
fi,o(z)  and  fz^i(z).  For  this  illustration,  the  number  of  video 
samples  integrated,  N,  is  chosen  to  be  nine  and  the  signal-to- 
noise  ratio  in  the  IF  range  bin  is  assumed  to  be  3  dB  when 
there  is  a  target  return.  The  choice  of  these  parameters  is 
in  part  to  facilitate  comparison  with  Marcum's  plots  [Fig. 4, 
Ref.  2]  which  assumed  infinite  video  bandwidth.  Accordingly, 
the  plots  in  Figure  3  for  y-O  are  found  to  be  in  agreement 
with  Marcum's  results  which  serves  as  a  validation  of  the 
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method  and  routines  developed  to  compute  the  probability 
density  functions.  Graphically,  the  effect  of  the  finite 
video  bandwidth  is  observed  to  be  two-fold.  Firstly,  the 
variance  of  tz,o(z)  is  reduced  vis-a-vis  the  density  function 
for  infinite  video  bandwidth.  Secondly,  the  variance  of 
is  also  reduced  but  this  is  accompanied  by  a  leftward 
shift  of  the  mode.  The  net  result  is  that  for  the  same 

probability  of  false  alarm,  the  detection  performance  has 
degraded  considerably.  This  is  the  consequence  of  the 
integration  of  additional  noise-only  samples  induced  by  the 
finite  bandwidth  of  the  video  amplifier. 

Figure  4  shows  the  receiver  operating  characteristics  for 
Y-0  and  Y"1*  The  number  of  video  samples  integrated  N  is 
equal  to  64.  For  a  given  signal-to-noise  ratio  and 
probability  of  false  alarm,  it  can  be  seen  that  the 
probability  of  detection  achievable  for  the  system  with  the 
finite  video  bandwidth  (Y"1)  is  poorer  than  the  system  with 
the  theoretical  reference  of  infinite  video  bandwidth. 
Conversely,  the  system  with  the  finite  video  bandwidth  will 
require  a  higher  signal-to-noise  ratio  to  have  the  same 
performance  as  the  system  with  infinite  video  bandwidth.  This 
is  illustrated  in  Figure  5.  For  a  probability  of  false  alarm 
and  probability  of  detection  equal  to  10“’  and  0.95 
respectively,  a  system  with  infinite  video  bandwidth  will 
require  an  input  signal-to-noise  ratio  of  0  dB.  To  achieve 
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Probability  of  False  Alarm 


the  same  performance,  the  system  with  a  video  bandwidth  that 
is  half  the  IF  bandwidth  will  require  an  input  signal-to-noise 
ratio  of  0.61  dB.  From  Eqn  (1),  this  implies  a  collapsing 
loss  of  0.61  dB  that  is  attributable  to  the  'insufficient' 
bandwidth  of  the  video  amplifier. 
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III. 


COLLAPSING  RATIO 


The  statistical  analysis  in  the  preceding  section  provides 
the  theoretical  background  and  insight  concerning  bandwidth 
effects  of  a  Gaussian-shaped  video  filter  on  the  performance 
of  the  square-law  detector.  The  analysis  has  also  shown  how 
sets  of  receiver  operating  characteristics  such  as  those  shown 
in  Figure  4  can  be  generated  using  a  digital  computer  to 
determine  the  input  signal-to-noise  required  to  attain  a 
desired  probability  of  detection  for  a  specified  probability 
of  false  alarm  and  a  given  IF  bandwidth  to  video  bandwidth 
ratio.  While  such  an  approach  is  feasible,  it  is  more 
convenient  in  practice  to  follow  Marcum's  approach,  using  the 
collapsing  ratio  p  defined  in  Eqn  (2)  to  determine  the 
required  input  signal-to-noise  ratio.  This  is  because  the  use 
of  the  collapsing  ratio  will  obviate  the  need  to  generate  sets 
of  receiver  operating  characteristics  for  each  value  of  y  • 
As  explained  in  Section  I,  only  the  receiver  operating 
characteristics  for  an  infinite  video  bandwidth  (Y“0)  are 
necessary  once  the  collapsing  ratio  is  known.  Furthermore, 
there  are  usually  other  sources  of  collapsing  loss  apart  from 
insufficient  video  bandwidth  in  practical  systems.  In  such 
cases,  the  formulation  of  the  collapsing  ratio  induced  by  each 
source  will  facilitate  computation  of  the  overall  collapsing 
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loss.  A  general  treatment  of  the  problem  of  multiple  sources 
of  collapsing  loss  can  be  found  in  Meyer  [Ref.  4]  and  Barton 
[Ref  5].  It  is  however  beyond  the  scope  of  this  report.  In 
this  section,  the  aim  is  to  formulate  the  collapsing  ratio 
associated  with  insufficient  video  bandwidth  for  the  Gaussian 
square -law  system  and  to  compare  the  result  with  Barton ^s 
formula  for  the  collapsing  ratio  in  Eqn  (5). 

From  Marcum's  definition  of  the  collapsing  ratio, 

_  _  M  +  N 

P  “  ~N~  (39) 

where  M  is  the  number  of  additional  noise-only  samples 
integrated  along  with  the  W  signal-plus-noise  samples. 

From  the  analysis  of  Section  II,  the  number  of  additional 
noise-only  samples  integrated  for  every  signal-plus-noise 
sample  integrated  in  the  Gaussian  square  law  system  is  equal 
to  Ny-1  where  Ny  is  given  by  Eqn  (26).  It  follows  that  if  the 
number  of  signal-plus-noise  samples  integrated  is  equal  to  W, 
then 


M  -  Nx(Ny-l)  (40) 


and 


N(Ny-l)  +  N 

Pg  -  - -  -  Nv  (41) 
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Hence,  the  collapsing  ratio  caused  by  insufficient  video 
bandwidth  in  the  Gaussian  square-law  system  is  equal  to  the 
equivalent  number  of  independent  IF  samples  averaged  by  the 
video  filter. 

The  subscript  'G'  has  been  used  to  denote  that  the 
collapsing  ratio  pg  given  in  Eqn  (41)  is  caused  by  insuffi¬ 
cient  video  bandwidth  in  a  Gaussian  square-law  system.  It  is 
also  used  to  distinguish  from  Barton's  formula  for  the 
collapsing  ratio: 


Pb 


Bjf  ^  2B^ 

2B„ 


1  +  Y 


(42) 


For  convenience,  pg  and  p^  shall  respectively  be  referred 

as  the  Gaussian  system  formula  and  Barton's  formula  in  the 
context  of  collapsing  ratios  caused  by  video  bandwidth.  As 
pointed  out  in  Section  I,  Barton's  formula  models  the  actual 
filters  with  ideal  rectangular  shaped  filters  having  band- 
widths  equal  to  the  noise-equivalent  bandwidths  of  the  actual 
filters.  A  comparison  of  the  two  formulae  will  provide  a 
relative  measure  of  the  accuracy  of  Barton's  formula  in 
modelling  systems  that  approximate  the  Gaussian  square-law 
system . 

Figure  6  plots  the  collapsing  ratio  formulae  pg  and  Pj,  as 
a  function  of  the  ratio  y  of  IF  bandwidth  to  twice  the  video 

bandwidth.  The  plot  obtained  by  applying  Barton's  formula  is 
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Ratio  of  IF  bandwidth  to  twice  video  bandwidth 

Figure  6.  Collapsing  ratio  due  to  insufficient  video 
bandwidth . 


indicated  in  the  figure.  As  eitpected  from  Eqn  (42),  it  is  a 
linear  plot  and  is  independent  of  input  signal-to-noise  ratio. 
All  the  other  curves  are  obtained  by  applying  the  Gaussian 
system  formula  for  input  signal-to-noise  ratios  of  0,  0.1,  l, 
and  2.  It  is  is  seen  that  is  a  function  of  the  input 

signal-to-noise  ratio  X.  For  the  same  y  /  the  collapsing 
ratio  given  by  p,,  decreases  as  the  signal-to-noise  ratio 

increases.  This  dependency  complicates  the  method  of  using 
the  collapsing  ratio  to  determine  the  required  signal-to-noise 
ratio  to  attain  a  specified  probability  of  detection  and 
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probability  of  false  alarm  since  the  collapsing  ratio  is 
itself  a  function  of  the  signal-to-noise  ratio.  For  very 
small  signal-to-noise  ratios  the  collapsing  ratio 

can  be  approximated  as 
Ps  -  Vl  +  2y' 

which  is  independent  of  the  input  signal-to-noise  ratio. 

The  collapsing  ratio  given  by  Eqn  (43)  can  be  considered 
as  the  worst  case  situation  for  any  given  y  in  the  Gaussian 
square-law  system.  Taking  pg-^/l+2Y^  as  the  reference,  it  is 
seen  that  Barton's  formula  overestimates  the  collapsing  ratio 
for  Y  less  than  two.  In  this  region,  the  largest  relative 
deviation  is  about  22.5%  in  the  vicinity  of  For  y>2, 

Barton's  formula  underestimates  the  collapsing  ratio.  For 
large  y ,  Ps/Pg  approaches  l/yf2  so  that  Barton's  formula 
underestimates  the  collapsing  ratio  by  a  factor  of  about 
for  very  large  y  • 

The  collapsing  loss  associated  with  each  collapsing  ratio 
can  be  determined  by  finding  the  incremental  increase  in 
signal-to-noise  ratio  required  with  respect  to  the  signal-to- 
noise  ratio  required  in  the  ideal  case  of  infinite  video 
bandwidth  for  the  same  detection  performance.  Marcum  has 
found  that  the  collapsing  loss  is  a  function  of  the  number  of 
signal-plus-noise  samples  integrated  but  varies  only  slightly 
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with  the  probability  of  detection  and  probability  of  false 
alarm  [Ref.  2].  His  results  for  the  collapsing  loss  as 
a  function  of  the  collapsing  ratio  y  and  with  W  as  a  parameter 
is  shown  in  Figure  7  for  P^-0.5  and  P^^-10‘“. 


Collapsing  rolio  p 


Figure  7.  Collapsing  loss  vs.  collapsing  ratio  for 
Pj,-0.5  and  Pf,-10'“.  (From  Ref.  4) 


Using  Figure  7  and  the  respective  equations  for  and  p^, 
the  collapsing  losses  predicted  by  the  application  of  Barton's 
formula  are  compared  against  those  deduced  using  the  Gaussian 
system  formula  for  selected  values  of  y  in  Table  I  and  Table 
II.  The  collapsing  loss  associated  with  the  use  of  Barton's 
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formula  is  denoted  in  the  table  by  L^.b  while  that  associated 
with  the  Gaussian  system  formula  is  denoted  by  Lc,,*  Using 
as  the  reference,  the  difference  between  the  two  collapsing 
loss  for  a  given  y  is  denoted  by  In  Table  I,  the  input 

signal-to-noise  ratio  X  is  assumed  to  be  much  less  than  one  so 
that  the  collapsing  ratio  can  be  taken  to  be  equal  to 
y/l  +  2y^  .  A  signal-to-noise  ratio  equal  to  one  is  assumed  in 

Table  II  to  provide  a  relative  bound  on  the  collapsing  loss 
for  signal-to-noise  ratios  between  zero  and  one. 

The  comparison  shows  that  for  small  signal-to-noise 
ratios,  the  collapsing  loss  predicted  by  Barton^s  collapsing 
ratio  is  within  0.3  dB  of  that  deduced  by  the  Gaussian  system 
formula.  If  the  signal-to-noise  ratio  is  near  or  equal  to 
one.  Barton's  formula  overestimates  the  collapsing  loss  by 
about  0.4  dB  for  video  bandwidths  around  half  the  IF 
bandwidth . 

Typically,  the  number  of  pulses  integrated  in  a  search 
radar  may  range  from  thirty  to  forty.  In  such  cases,  the 
required  single-pulse  input  signal-to-noise  ratio  is  close  to 
unity  for  reliable  detection  performance  [Ref.  5].  From 
Table  II,  it  follows  that  the  collapsing  loss  will  be 
overestimated  by  about  0.4  dB  if  Barton's  formula  is  used  in 
a  Gaussian  square-law  system  with  a  video  bandwidth  about  half 
the  IF  bandwidth. 
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TABLE  I.  COMPARISON  OF  COLLAPSING  LOSS. 
BARTON'S  FORMULA  VS  GAUSSIAN  SYSTEM  FORMULA. 


Input  SNR 

X  «  1 

Collapsing  loss,  L.  dB 

N  •  10 

N  ■  1000 

B 

Pm 

Po 

D 

'■c.g 

B 

B 

1.0 

m 

1.7 

D 

m 

m 

1.1 

0.8 

B 

m 

H 

B 

1.6 

1.8 

m 

2.5 

2.8 

B 

m 

6.0 

7.1 

2.2 

B 

m 

3.3 

3.6 

B 

TABLE  II.  COMPARISON  OF  COLLAPSING  LOSS. 
BARTON'S  FORMULA  VS  GAUSSIAN  SYSTEM  FORMULA. 


Input  SNR 

X  -  1 

Collapsing  loss,  L.  dB 

V 

O 

■ 

z 

N  •  1000 

B 

B 

B 

B 

'■c.g 

B 

1.0 

B 

1.5 

D 

B 

B 

1.1 

B 

B 

B 

2.6 

1.2 

1.0 

B 

1.8 

B 

B 

B 

6.5 

2.2 

2.3 

-0.1 

3.3 

B 

-0.1 
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IV.  CX)NCLUSIONS 


The  effects  of  Gaussian-shaped  IF  and  video  filters  on  the 
performance  of  a  square-law  detector  in  pulsed  radar 
applications  employing  post-detection  integration  have  been 
analyzed  in  this  paper. 

Specifically,  a  formula  for  determining  the  number  of 
additional  noise-only  samples  integrated  due  to  the  finite 
bandwidth  of  the  video  filter  has  been  derived.  Emerson's 
formulation  for  the  cumulants  of  the  voltage  probability 
density  function  at  the  output  of  the  video  filter  has  been 
used  in  this  derivation. 

The  resultant  cumulants  of  the  probability  density 
function  at  the  output  of  the  post-detection  integrator  are 
obtained.  Based  on  these  cumulants  and  Edgeworth's  asymptotic 
expansion  of  the  Gram-Charlier  series  to  represent  the  density 
functions,  MATLAB  routines  have  been  written  to  plot  the 
receiver  operating  characteristics  for  various  ratios  of  IF 
bandwidth  to  video  bandwidth.  The  degradation  in  performance 
or  collapsing  loss  can  be  observed  and  quantified  from  these 
operating  characteristics  as  the  video  bandwidth  is  reduced 
relative  to  the  IF  video  bandwidth. 

Taking  into  account  the  Gaussian  shape  of  the  filters, 
the  collapsing  ratio  is  found  to  be  equal  to  the  equivalent 
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number  of  independent  IF  samples  averaged  by  the  video 
filter.  This  is  given  by  Eqn  (25)  which  is  a  quadratic 
equation  in  with  signal-to-noise  ratio  and  IF  bandwidth  to 
video  bandwidth  ratio  as  parameters.  The  collapsing  ratios 
given  by  Eqn  (25)  and  the  associated  collapsing  loss  are 
compared  against  those  predicted  by  Barton^ s  formula  in  Table 
I  and  II.  The  comparison  shows  that  for  typical  signal-to- 
noise  ratios  close  to  0  dB  and  video  bandwidths  about  half  the 
IF  bandwidth.  Barton ^s  formula  for  the  collapsing  ratio 
results  in  an  overestimation  of  the  collapsing  loss  by  about 
0.4  dB. 
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APPENDIX  A 


MATHEMATICAL  DERIVATION  OF  THE  CUMULANTS 
This  appendix  follows  the  method  of  Emerson  [Ref.  7]  in 
deriving  a  closed  form  expression  for  the  cumulants  of  the 
probability  density  function  f^(Y)  at  the  output  of  the 
Gaussian  video  filter  shown  in  Figure  A-1. 


Figure  A-1.  Block  diagram  of  the  receiver. 


The  input  to  the  system  is  a  Gaussian  random  process  given 
by: 

( t)  -  s(  t)  +  n(  t)  (A-1) 

where 

S(t)  -  coso^t  (A-2) 

and  n(t)  is  additive  white  Gaussian  noise  (AWGN)  with  a  power 
spectral  density  equal  to  r\/2. 


The  impulse  response  functions  of  the  Gaussian-shaped  IF 


filter  and  video  filters  are  respectively  defined  by: 

JL  f.2 

h^f(t)  -  2  (2jcp2)  2  [  _  ^27iP)  ]  cos2nf^t  (A-3) 

and 


— 

-  (2tcv2)  2  exp  [  -  (2nv)  2__  ]  (A-4) 

where  P  and  v  are  respectively  the  root-mean-square  (rms) 
bandwidths  of  the  IF  and  video  filters. 

The  starting  point  in  this  derivation  is  the  following 
general  expression  of  the  cumulants  Kr  obtained  by  Emerson  for 
an  arbitrary  input  signal  contaminated  by  AWGN  into  a  square- 
law  detector  with  IF  and  video  filters  of  arbitrary  frequency 
response  characteristics: 

K,-  ’  EM  ^  (A-5) 

^  i-O  i-0 

where 

«o 

s^{t)  -  J"  s(t-x)  ^^(x)  dx  (A-6) 
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Xj  and  are  respectively  the  eigenvalues  and  eigen¬ 


functions  defined  by  the  integral  equation: 


Xj<|)^(x)  -  J sr(x,y)  4>^(y)  dy 


(A-7) 


The  function  g(u,v)  is  called  the  system  kernel  and  it  is 
completely  defined  by  the  impulse  response  functions  of  the 
filters.  In  particular,  it  is  defined  by  the  following 
equation: 

giu.v)  -  h^(z)  h^f(v-z)  dz  (A-8) 

<-00 

which  can  be  expanded  into  the  uniformly  convergent  bilinear 
Mercer  series,  i.e., 

00 

g{u,v)  -  52  (A-9) 

i-O 

To  determine  the  cumulants  for  the  Gaussian  system  shown 
in  Figure  A-1,  the  method  of  solution  is  to  evaluate  the 
system  kernel  using  Eqn  (A-8)  and  expressing  it  in  a  form  that 
matches  the  Mercer  series.  The  eigenvalues  and  eigenfunctions 
that  satisfy  the  integral  equation  of  Eqn  (A-7)  can  then  be 
identified.  The  cumulants  can  then  be  determined  by  using 
these  eigenvalues  and  eigenfunctions  to  evaluate  Eqn  (A-5). 
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r 


Thus  using  Eqn  (A-8),  the  system  kernel  for  the  Gaussian 
system  shown  in  Figure  A-1  is  given  by 

«• 

g{u,v)  -  j"  2yzirpexp[-(2TEfi)^  ]  cos<»>o(u-z) 

2 

X  yjiTvexp  [  -  (2iiv)  2 -^  ] 

X  2v'2irpexp  [  -  (27tp)  2 ]  cosci>j,(  v-z) 

X  dz  (A-lO) 


Now, 


COS«<,(u-z)  COS<iig(V-z) 


(A-11) 


-  -^COS<i)g(u-V)  + 


[  sin<«)j,(u+v)  sin2(i>oZ  +  sinci)p(u+v)  sin2w^z  j 


Substituting  Eqn  (A-11)  into  Eqn  (A-10), 
g{u,v)  -  4xp^>/27r  [gr^  +  g^]  (A-12) 
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where 


-  cos<i>o(u-v)  J  exp  [  -  {  J?JLPL-  [  (u-z)^  +  (v-z)^] 


+  }  1 
2  J 


92  ~  j  cos(i>^(u+v-2z)  exp  [  -  {  [  {u-z)^+  (v-z)^] 


j  dz 


(A-13) 


2  -I 


(A-14) 


The  integral  in  Eqn  (A-13)  can  be  evaluated  by  completing 
the  square  in  the  exponent  and  casting  the  integrand  into  the 
form  of  the  standard  integral  identity:  [Ref.  12] 


(A-15) 


ca  _ 

j  exp[-a^  (x-m)^]  dx  - 


Thus  completing  the  square. 


-  COSO)  (u-v)  exp  f  {u^  +  v^)  -  }  | 

^  ^  2  1+2y^  1+2y^  ^ 


oo 

f  exp  {  (1+2y^)  Iz - ^ — r  (u+v)  }  dz 

J  1+2y2 


(A-16) 
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where 


1  liL 

V  2B^ 


(A-17) 


The  symbols  and  B,  in  the  above  equation  are 
respectively  the  noise  equivalent  bandwidths  of  the  Gaussian 
IF  and  video  filters. 


Making  use  of  Eqn  (A-15)  with  a 


y2 


(1+2y2)  , 


COS<i>q(U-V)  exp  [  (u2  +  v2) 

o  ••  4  1+2y^ 


v/2irv  (1+2y^)  ^ 


1+2y^  ' 
(A-18) 


Now,  consider  the  second  term  ga  of  the  system  kernel 
given  in  Eqn  (A-12).  By  completing  the  square  in  the  exponent 
as  in  the  evaluation  of  g^,  it  can  be  similarly  shown  that 


3,  . 

4  1+2Y^  1+2y^ 


(A-19) 


X  j  cosu)g{u+v~2z)  exp  [-a^(z-n?)^]  dz 


where 


a  - 


m  - 


(1  +  2y2)  2 
Y 


— (u+v) 
1  +  2y^ 


(A-20) 


(A-21) 
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Now, 


j  cos(iig(u+v-2z)  exp  [-a^  iz-in)^]  dz 

(A-22) 

-  J  cosUj,(i2+v)  cos2&)o2  exp  [-a^  (z-jn)  ^  ]  dz 

«a 

+  J  sin<«)o(u+v)  sin2tt>j,z  exp  [ -a^  (z-m)  2]  dz 

—  09 

Let  X  -  (z-m)  .  Then, 

«a 

J  coso)q(u+v)  cos2(«>^z  exp  [ -a^  (z-m)  dz 

(A-23) 

flO 

-  cos<i)o(u+v)  cos2Wom  J  exp[-a^x^]  dx 

-«o 

flO 

cosfo^{u+v)  sin2(t>jn  j  sin2(i>^  exp [ -a^x^ ] 

— «o 

dx 

From  Eqn  (679)  of  Ref.  12, 

/*  cos2<i)_x  exp  [ -a^x^]  dx  -  exp  [  -  ] 

J  3  " 

-aa 

(A-24) 

«• 

J  sin2<i>^  exp  [ -a^x^  ]  dx  -  0 

-•a 

(A-25) 
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Hence , 


J  cos<i>^(u+v)  cos2(o^z  exp[-a2(z-m)2]  dz 


i —  f  Ci)^ 

-  cos(i)o(u+v)  cos2<i)^  -5^  exp - ^ 

a  ^  a^ 


(A-26) 


The  second  integral  in  Eqn  (A-22)  can  be  similarly 
evaluated  to  give  the  following  result: 

m 

J  3in<i)g(u+v)  sin2<i)oZ  exp  [ -a ^  ( z-;n)  ^ ]  dz  (A-27) 


2 

sin«^(u+v)  sin2<o^  exp  I 


Using  the  results  of  Ern  v\-22),  Eqn  (A-26),  and 
Eqn  (A-27)  to  evaluate  Eqn  (A-19),  it  can  be  shown  that 


9^,  .  exp[--«^( 

4  1+2y2  1+2y^ 


(A-28) 


X  exp  [ 


cos 


2Tc2p2  (2+2y^) 


{  (i)  (u+v)  ^  } 

_ l-t-2Y^ 

y/2n  V  (1+2y^)  ^ 


Comparing  Eqn  (A-28)  with  Eqn  (A-18),  it  is  obvious  that 
if  o>o  >  2xP,  then  will  be  small  and  negligible  compared  to 
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g^.  Hence,  neglecting  the  contribution  from  g^,  the  system 
kernel  given  in  Eqn  (A--12)  becomes 


giu,  v) 


47CP^COS  WqCu-v) 

v'1+2y^ 


(A-29) 


X 


exp  ^ 


r  _  (2mP)^ 


1+2y' 


(u^+v^) 


2y^u V}  ] 
1+2y^  ^ 


Using  the  following  substitutions. 


(1+2y^)  ^ 

u  -  ^  —  X 

2x3 

„  .  UliY!) 


1 

2\  4 


2x3 

y/l+2y^  -1 
yjl+2y^  +1 


(A-30) 


(A-31) 


(A-32) 


Eqn  (29)  can  be  transformed  into  the  following  form  suitable 
for  the  application  of  Mehler's  formula: 

g(u,v)  -  [  cos <o-u  cos u_v  +  sinw^u  sinu.v  ]  (A-33) 


X  exp  [  -  ]  exp  [  -  -  2t:xy 


Using  Mehler's  formula. 


[  tMx=^y2)  -2txy  ]  .  yjrp-  f  ^ 


(A-34) 
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where  H^(x)  is  the  Hermite  polynomial  defined  by 


-  (-1)-^  exp[x^]  expL-x^] 

dx^ 


(A-35) 


From  Eqn  (A-30)  and  Eqn  (A-31), 

(A-36) 

(A-37) 

Substituting  the  above  equations  into  Eqn  (A-33), 


X  - 


2nfiu 


y  - 


(1+2y^)  * 
2ii3v 

(1+27^)  4 


g(u,v)  -  g^{u,v)  +  g^{u,v) 


(A-38) 


Where 


ZpyiT  f  ^/l-^^2Y^  -1  Ti 

1+Vi+2y^  }/l+2y^  +1 


X 


X 


X 


y4xP  exp  [  -  ]  cos  WpU 

_ \/l+2y^  ^ _ 


[  2^i!  V^(1+2y^)  ^  ^ 


y/Urp"  exp  I  -  ■  —  P-  -  ——  j  cos 
_ yjl^2y^  2 _ 

±  j, 

[  2^i!  ^/W  (1+2y^)  *  ]  ^ 


±  ]  «,  [ 

(1+2y^)  *  (1+2y^)  * 


(A-39) 
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g^^U,v)  -  ]i 

1+Vi+2y^  yjl+2y^  +1 


(A-40) 


V4itp  exp  [  j  cosw^u 

_ ^1+2y^  2 _ 

^  _i 

[  2^i!  yir(l+2Y^)  *  3  2 


>/4nft  exp  ^ P ^ I  sinWgV 
_ )/l+2y^  ^ _ 

j,  _i 

[  2^i!  ^  (1+2y^)  ^  3  2 


X 


2xP  u 

J. 

(1+2y^)  ^ 


1  „  r  27tp  V 

L - ^ 

(1+2y^)  ^ 


By  comparing  Eqn  (A-38)  against  the  Mercer  series 
expansion  of  the  system  kernel 


g^{u,v)  -  T 


(A-41) 


the  eigenvalues  and  eigenfunctions  can  be  deduced  to  be 


2Pv^  f  ^l+2y^  -1  ii 

1+  Vi+2y^  /T+Ip’+i 


(A-42) 


^  exp  [  - ]  H,  [  ^  ]  cos»,u 


1  1 


[  2^i!  ^(1+2y2)  *  3  2 


(A-43) 
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similarly,  it  can  be  shown  that 

^s,i  “  ^c,  i 


(A-44) 


nr^mr _ T  (2710)^  1  „  f  27i3u  1  _ _ 

v4Wp  exp  I  —  r —  __  j  1^  - ^ — -  j  sinWgU 

_ (1.2Y^)  ‘ _ 


J,  _1 

[  2^i!  ypK  {1+2y2)  ]  2 


(A-45) 


Hence,  from  Eqn  (A-42)  and  Eqn  (A~44), 


Z  -  t- 


(A-46) 


i-0 


i-0 


(2Pv/^)  ^  ^  f  -1  ^  ri 

(1  +V^1+2y^  )  ^  \/i+2y^  +1 


Now, 


Y  - —  for  I  X  I  <  1 

h  1-^ 


Hence , 


i-0  (Vi+2y^+1)'"-  (v'1+2y2-1)^ 


i-O 


(A-47) 


E  -  E  *  t 


i^o 


1-0 


i-0 


2  (2pv^)^ 


(>/i+2y^+1)^-  Wl+2y^-l)^ 


(A-48) 
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From  Eqn  (A-12)  and  Eqn  (A-45), 


-  yf2A  j  cosa)^( t-u)  du 


v/2~A  (4itp)  ^  (A-49) 

[  2^ily/ir  (l+2y^)  *  ]  2 

aa 

X  j  exp(-a^-^)  H^iau)  cosu^u  cosa)„( t-u)  du 


where 


2irP 

_i 

(1+2y^)  ^ 


(A-50) 


1  1 

cos  WqU  COSWo(  t-u)  -  -i-COSWpt  +  -icos  {20>pU+ t) 

(A-51) 

For  Ci)o>2iiP,  the  residue  from  the  high  frequency 

component  in  Eqn  (A-51)  can  be  neglected  in  the  integration  of 
Eqn  (A-49).  Hence, 


vTA(4nP)^  COSto^t  (A-52) 

2.  1  2 

[  2^ily/n  (l+2y^)  «  ]  2 

«a 

X  J  exp(-o^-^)  H^iau)  du 
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The  following  integral  identities  derived  in  Ref.  13  are 


useful  in  evaluating  (A-52): 


j  exp  (  -x^ )  ( ^P^x)  dx  *  0 


(A-53) 


0* 

J  exp(-x2)  dx  - 


(A-54) 


Let  “  x^ 

2 


(A-55) 


Then  au  -  ^x  and  du  -  — dx. 

a 


Using  the  above  substitution,  Eqn  {A-53)  and  Eqn  (A-54) 


J  exp(-a^^)  H^iau)  du 


/  exp(-x2)  H,{y/2x)  dx 


^  y/2n  (1+2y^)  ^  i! 

I  2xp  , 

2 


for  i  odd 


for  i  even 


(A-56) 
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Hence,  using  the  result  of  Eqn  (A-56)  to  evaluate  Eqn  (A-52) 


r  (  \l2A^  j  2 

I  2pv^ 


(1+2y^)  ® 


(i!)  ^ 


COSlDgt 


for  1  even. 


s 


C,  1 


(t)  -  <1 


(A-57) 


0  for  i  odd. 


A  similar  result  is  obtained  for  s.^ift;  except  that  the 
cosUot  term  in  Eqn  (A-57)  is  replaced  by  a  sinu^t  term. 

Using  Eqn  (A-40)  and  Eqn  (A-55),  it  can  be  shown  that 

V  (C)  -  (2Pv^)_^  (1+2y^)  ^  cos^Q^t 

-  '  (l+yi+2Y2 

^  r  V^H-2Y^  -1  I2ri  (2i)  ! 

(A-58) 

Let  [  -  X  (A-59) 

/r+2Y^+i 


From  Ref.  12, 


j,i  (2i)  1 


—  (  — +  I)x2 

1  +  Ax  ^  2  2  ^ 

2  2  ; 


-  (1-x)  2 


(A-60) 
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Hence,  using  this  result  to  evaluate  Eqn  (A-58)  gives 


t,  sj.l(t) 


2A^(2Pv/tc  {l+2y^) 


2-0 


[  (V1+2y^+1)2^-  (Vi+2y^-1)^^  ]  2 


(A-61) 


Similarly,  it  can  be  shown  that 


2A^  (2^y/W)^-^  (1+2y^)  * 

_i 

[  (V1+2y^+1)^^-  (V1+2y^-1)^^  ]  2 


sin^w^t 


Hence , 


(A-62) 


«• 


E  iK.i)"  4i(t)  ^  E 


_1 

2A^  (2P>/7r)  (1-*-2y^)  ^ 

_i 

[  (V1  +  2y^  +1)^^  -  {/i+2y^  -1)^^  ]  ^ 


(A-63) 
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From  Eqn  (A-5)  and  the  results  of  Eqn  (A-48)  and  Eqn  (A-63), 


K,  -  sl(t) 

2  i-0  i-o 

_  (2TiPvD^(r-i)  !  ^ 

(v/i+2y^  +1)^  -  (V'1+2y^  -!)"■ 


r 

L 


1  +  r- 


N 


(^/1-«-2y^ -«-l)  (v/1-«-2y^  -1)^ 

[  (v^1+2y^+1)^-  (/1+2y^-1)^  ]  ^ 


X 


(1+2y=) 


(A-64) 


Now,  it  can  be  shown  that  the  IF  noise  bandwidth  is 
equal  to  P>/Tr.  It  follows  that  the  total  average  power  into 
the  detector  is  given  by 

-  TiP\/ir  (A-65) 


Hence , 


_ (2o^)^(r-l)  ! _ 


(A-66) 


r 

i 


1 


+  r 


A‘ 


N 


(Vi-*-2y^-*-i)^  -  (^Jl+2y^-l)^ 
(>/i  +  2y^+1)''  +  (v^iT2y^-1)^ 


(1+2y") 
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This  equation  which  expresses  the  cumulants  in  terms 
of  the  IF  bandwidth  to  video  bandwidth  ratio  y  and  the  input 
signal-to-noise  ratio  A^/a^  concludes  the  mathematical 
derivation  of  this  appendix. 
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APPENDIX  B 


HATLAB  PROGRAMS  FOR  PLOTTING 
THE  RECEIVER  OPERATING  CHARACTERISTICS 

This  appendix  is  a  listing  of  the  Matlab  programs  written 
as  part  of  this  thesis  to  plot  the  probability  density 
functions  and  operating  characteristics  of  the  Gaussian  square 
law  system  described  in  Section  II.  The  programs  are  written 
as  MATLAB  M-Files.  The  following  are  respectively  the 
filenames  of  the  programs  for  plotting  the  probability  density 
functions  and  operating  characteristics; 

1.  SQ2PDF,M:  This  plots  the  probability  density  functions 
based  on  Eqn  (29)  which  is  an  Edgeworth's  asymptotic 
expansion  of  the  Gram-Charlier  series. 

2.  SQ2ROC.MX  This  plots  the  receiver  operating  characteris¬ 
tics  based  on  Eqn  (37)  which  is  derived  by  integrating 
Eqn  (29)  from  the  threshold  voltage  to  infinity. 

These  are  script  files  which  are  run  by  invoking  their 
filenames  in  the  MATLAB  operating  environment.  Both  programs 
will  prompt  the  user  to  enter  the  necessary  input  parameters 
prior  to  execution.  In  addition  five  function  type  M-Files 
are  specifically  created  to  facilitate  programming  and 
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execution  of  these  script  files.  The  filenames  of  these  five 
functions  and  their  purpose  are: 

1.  GCUMU.Mz  This  function  computes  the  cumulants  of  the 
video  filter  output  probability  density  funtion  based  on 
Eqn  (10). 

2.  GDFZ,M:  This  function  computes  the  probability  density 
for  given  normalized  values  of  the  Gaussian  square- law 
system  output  based  on  Eqn  (29). 

3.  UPFZ.Mz  This  function  computes  the  probability  of  false 
alarm  or  probability  of  detection  of  the  Gaussian  square- 
law  system  for  a  specified  threshold  voltaqe.  The 
computation  is  based  on  Eqn  (37)  and  Eqn  (38). 

4.  GCR.Mz  This  function  computes  the  collapsing  ratio 
caused  by  the  finite  video  bandwidth.  It  is  based  on  Eqn 
(26). 

5.  HMT.Mz  This  function  computes  the  mth  order  Hermite 
polynomial  for  values  of  m  less  than  or  equal  to  12. 


The  MATLAB  source  programs  of  the  two  script  files  and  the 
five  function  files  are  provided  in  the  remaining  pages  of 
this  appendix. 
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1.  SQ2PDF.M  PROGRAM 


%  ========================================================== 

%  This  program  computes  and  plots  the  pdf  of  a  square  law 
%  detector  with  Gaussian  IF  and  video  filters.  It  will 
%  prompt  the  user  to  enter  the  following  parameters: 

%  a.  Range  and  resolution  of  the  x-axis  for  the  plot. 

%  b.  SNR  =  input  signal-to-noise  ratio.  (Not  in  dB) 

%  c.  GAMMA  =  IF  bandwidth  to  twice  video  bandwidth  ratio. 

%  d.  N  =  Number  of  pulses  to  be  integrated. 

%  In  addition  to  the  pdf  plot  for  the  specified  GAMMA,  the 
%  pdf  plot  for  an  infinite  video  bandwidth  will  automatically 
%  be  plotted  for  comparison  purposes. 

% - 

disp( 'The  output  pdf  will  be  plotted  as  a  function  of  the' 
disp( 'normalized  voltage,  z.') 

disp( 'Enter  range  and  resolution  of  z  to  be  plotted') 
range  =  input( ' [zmin,zmax,zres]  :  ') 

N  =  input ( 'Number  of  pulses  to  be  integrated  =  ') 

SIGMA  =  input ( Input  noise  rms  voltage  =  ' ) 

SNR  =[00] 

GAMMA  =[00] 

ORDER  =  4 

SNR(2)  =  input ('Enter  input  SNR:  ') 

GAMMA(2)  =  input( 'Enter  value  of  GAMMA:  ') 
disp( 'Computing  ...') 

KYO  =  [] 

KYI  =  [] 

Z  =  range(l) :range( 3) :range(2) 

GAMMAO  =  GAMMA(l) 
for  j  =  1:2 
K  =  [] 

SNRj  =  SNR(j) 
for  r  =  1: ORDER 

K(r)  =  gcumu(r, SIGMA, SNRj, GAMMAO) 
end 

if  j  ==  1 

fzinfo  =  gdpz ( ORDER, N,Z, SIGMA, SNRj, GAMMAO) 
else 

fZinfl  =  gdpz ( ORDER, N,Z, SIGMA, SNRj, GAMMAO) 
end 
end 

GAMMAV  =  GAMMA(2) 

SNRj  =0.0 

for  r  =  l:ORDER+2 

KYO(r)  =  gcumu (r, SIGMA, SNRj, GAMMAV ) 
end 

fZvO  =  gdpz ( ORDER, N,Z, SIGMA, SNRj, GAMMAV, KYO) 
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SNRj  =  SNR(2) 

NPRIME  =  gcr(SNRj,GAMMAv) 

SNRjbar  =  SNRj /NPRIME 
for  r  =  l:Order+2 

KYl(r)  =  gcuittu(r, SIGMA, SNR jbar,GAMMAv) 
end 

fZvl  *  gdpz ( ORDER, N,Z, SIGMA, SNR j,GAMMAv, KYI) 

plot(Z,fZinfO, ,Z,fZvO, '  — ',Z,fZinfl, ,Z,fZvl, '  — M 
grid 

xlabel( 'Normalized  Output  Voltage,  z') 
ylabel( 'Probability  Density,  f(z)') 


2.  SQ2ROC.M  PROGRAM 


%  This  program  computes  and  plots  the  ROC  data  points  for 
%  the  Gaussian  square  law  system.  It  is  programmed  to  provide 
%  accuracy  up  to  the  4th  grouping  of  Edgeworth's  series.  To 
%  obtain  any  other  order  of  accuracy,  the  local  variable 
%  ORDER  which  is  set  to  4  and  the  function  file  gdfz.m  must 
%  be  modified  accordingly. 

% - 

disp( 'Choices: ' ) 

disp('l.  Computing  Pfa  and  Pd  given  a  normalized  threshold.') 
disp('2.  Pd  vs  Pfa  curves.') 
disp( '  ) 

CHOICE  =  input ( 'Select  one  of  the  above:  ') 
if  CHOICE  ==  1 
ORDER  =  4 
MORE  =  1 
while  MORE  ==  1 

N  =  input ( 'Number  of  samples  integrated  =  ') 

SNR  =  input ('Input  SNR  =  ') 

SIGMA  =  input ('Input  rms  noise  voltage  =  ') 

GAMMA  =  input( 'GAMMA  =  ') 

VT  =  input ( 'Normalized  threshold  voltage,  VT  =  ') 

disp( 'Computing  ...  ' ) 
for  r  =  l:0RDER+2 

KYO(r)  =  gcumu(r,SIGMA,o, GAMMA) 
end 

NPRIME  -  gcr( SNR, GAMMA) 

SNRjbar  *  SNRj/NPRIME 
for  r  =  l:ORDER+2 

KYl(r)  =  gcumu(r, SIGMA, SNRjbar, GAMMA) 
end 


56 


Pfa  =  upfz(ORDER, N,VT, SIGMA, KY0,1) 

Pd  =  upf Z( ORDER, N,VT, SIGMA, KYI, 1) 

disp('Do  you  wish  to  compute  Pd  vs  Pfa  for  another  VT?') 
MORE  =  input Type  1  if  YES,  0  if  NO  :  ') 

end 
else 

N  =  input ( 'Number  of  samples  integrated  =  ') 

GAMMA  =  input( 'GAMMA  =  ') 

disp( 'Enter  up  to  a  maximum  of  3  SNR  values  as  elements  of 
[]') 

SNR  =  input ('SNR  []  =  ') 

SIGMA  =  input ('Input  rms  noise  voltage  =  ') 

disp( 'Enter  range  &  resolution  of  z  to  be  used  in 
computing ' ) 

disp( 'the  ROC.') 

range  =  input( ' [ zmin, zmax, zres] ' ) 

Z  =  range ( 1 ) : range ( 3 ) : range ( 2 ) 
for  r  =  l:0RDER+2 

KYO(r)  =  gcumu(r, SIGMA, 0, GAMMA) 
end 

NPRIME  =  gcr( SNR, GAMMA) 

SNRjbar  =  SNR  ./NPRIME 
for  r  =  l:0RDER+2 

KYl(r)  =  gcumu(r, SIGMA, SNRjbar, GAMMA) 
end 

Pfa(j,:)  =  upf Z( ORDER, N,Z, SIGMA, KY0,0) 

Pd(j,:)  =  upfz(ORDER,N,Z, SIGMA, KYI, 0) 
if  length (SNR)  ==  1 
semilogx(Pfa,Pd, '-' ) 
elseif  length (SNR)  ==  2 

semilogx(Pfa( 1 , : ) ,Pd( 1 , : ) ,Pfa( 2 , : ) ,Pd( 2 , : ) ) 
elseif  length (SNR)  ==  3 

semilogx(Pfa(l, : ) ,Pd( 1 , : ) ,Pfa( 2 , : ) ,Pd( 2 , : ) ,Pfa(3, : ) ,Pd(3, : ) ) 
end 
grid 

xlabel( 'Probability  of  False  Alarm') 
ylabel( 'Probability  of  Detection') 
end 
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3.  GCUMU.N  FUNCTION 

function  K  =  gcuinu(r, SIGMA, SNR, GAMMA) 

%  This  function  computes  the  rth  cumulant  of  the  output  pdf 
%  of  a  square  law  detector  with  Gaussian  IF  and  video  filters. 
%  GAMMA  is  the  ratio  of  the  IF  bandwidth  to  twice  the  video 
%  bandwidth.  SIGMA  is  the  input  rms  noise. 

%  Reference:  Thesis  Eqn  (10). 

% - 

G12  =  1  +  2*GAMMA^2 
G12PR  =  (sqrt(G12)  +  l)^r 
G12MR  =  (sqrt(G12)  -  l)^r 
BNl  =  (2  *  SIGMA^2)^r 
RMl  =  r  -  1 
BN2  =  fac(RMl) 

BDl  =  G12PR  -  G12MR 
B  =  (BNl  *  BN2)/BD1 
CNl  =  BDl 

CDl  =  G12PR  +  G12MR 
C  =  sqrt(CNl/CDl) 

K  =  B  *  (1  +  r*SNR*G12*0.25) 


4.  GDFZ.M  FUNCTION 

function  fZ  =  gdf 2( ORDER, N,Z, SIGMA, SNR, GAMMA, K) 

%  This  function  computes  the  probability  density  for  specified 

%  values  of  the  Gaussian  square  law  system  normalized  output, 

Z. 

%  K  is  an  array  comprising  the  cumulants  of  the  pdf. 

%  Reference:  Thesis  Eqn  (29). 

% - 

C  =  [] 

for  j  =  1: ORDER 
SUMn  =0.0 
SUMd  =0.0 
for  m  =  1:N 

SUMn  =  SUMn  +  K(m,j+2)/N 
SUMd  =  SUMd  +  K(m,j+2)/N 
C(j)  =  SUMn/(N^( j/2)  *  SUMd-(( j+2)/2)) 
end 

fz  =  zeros ( 1, length( Z) ) 

SUM  =  zeros(ORDER,length(Z) ) 

X  =  (SIGMA*Z  -  N*K(l))/sqrt(N*K(2)) 

CONST  =  SIGMA/sqrt(N*K(2)) 
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PSIO  =  exp(-(Z  . ^2)/2)/sqrt(2*pi) 
for  j  =  1: ORDER 
if  j  ==  1 

SUM(1,;)  =  CONST  *  PSIO  .*(  1  +  hmt(Z,3)  *  C(l)/6  ) 
elself  j  ==  2 

SUM(2,:)  =  CONST  *  PSIO  .*(  hmt(Z,4)  *  C(2)/24  +  hint(Z,6) 

*  (C(l))"2/72  ) 

elseif  j  ==  3 

SUM(3,:)  =  CONST  *  PSIO  .*(  hmt( Z , 5) *C( 3 )/12  +  hmt(Z,7)* 
C(1)*C(2)/144  +  hint(Z,9)*(C(l))''3/1296  ) 

elseif  j  ==  4 

SUM(4,:)  =  CONST  *  PSIO  .*(  hmt ( Z , 6 ) *C( 4 ) /720  +  hrot(Z,8)* 
(  ( C( 2) ) "2/1152  +  C(l)*C(3)/720  )  +  hmt (Z, 10)* 
(C(l))"2  *  C(2)/1728  + 

hmt(Z,12)*( (C(l) ) "4/31104) 
end 

f z  =  f z  +  SUM( j , : ) 
end 


5.  UPFZ.M  FUNCTION 

function  UPTE  =  upf Z( ORDER, N,Z, SIGMA, ERROR) 

%  ========================================================== 

%  This  function  computes  the  Pfa  or  Pd  of  a  Gaussian  square 
%  law  detector  for  a  specified  threshold  normalized  by  the  rms 

%  noise  voltage  at  the  input  to  the  detector.  If  the  argument 
%  ERROR  is  set  to  1,  the  function  will  return  a  two  element 
%  row  vector.  The  second  element  indicates  the  order  of 
%  accuracy  of  the  Pfa  or  Pd  value  given  by  the  first  element. 

%  If  ERROR  is  set  to  0,  then  the  function  will  return  the 

%  the  value  of  Pfa  or  Pd  without  any  indication  of  the  order 
%  of  accuracy. 

%  Reference:  Thesis  Eqn  (37)  and  Eqn  (38). 

% - 

C  =  [] 

for  j  =  1: ORDER 
SUMn  =  O.o 
SUMd  =0.0 
for  m  =  1:N 

SUMn  =  SUMn  +  K(m,j+2)/N 

SUMd  =  SUMd  +  K(m,2)/N 

end 

C(j)  =  SUMn/(N"( j/2)  *  SUMd"( ( j+2)/2) ) 
end 

ZZ  =  (SIGMA*Z  -  M*K(1) )/sqrt(M*K(2) ) 

PSIO  =  exp(-ZZ  . "2)/sqrt(2*pi) 

ZM  =  ZZ/sqrt(2) 

SUMl  =  0.5*erf (-INF, ZM, 'high' ) 
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TERMl  =  hint(ZZ,2)*C(l)/6 

TERM2  =  hint(ZZ,3)*C(2)/24  +  hmt  (  ZZ  ,  5  )  *C(  1 )  ^2/72 
TERM3  *  hmt(ZZ,4)*C(3)/120  +  hrot( ZZ , 6 ) *C( 1) *C( 2 )/144  + 
hint(ZZ,  8)  *0(1)^3/1296 

TERM4  =  hmt(ZZ,5)*C(4)/120  +  hmt ( ZZ , 7 ) *C( 2 ) ^2/1152  + 

hint(ZZ,7)*C(l)*C(3)/720  +  hmt  (  ZZ ,  9)  *C(  1 )  ^2*C(  2  )/1728 
hiat(ZZ,  11) *0(1)  ^4/31104 
SXJM2  =  PSIO  .*  (TERMl  +  TERM2  +  TERM3  +  TERM4) 
if  ERROR  ==  1 

UPTE(l)  =  1  -  (SUMl  -  SUM2) 

UPTE(2)  PSIO  .*TERM4 
else 

UPTE  =  1  -  (SUMl  -  SUM2) 
end 

%  ========================================================== 


6.  GCR.N  FUNCTION 
function  GRHO  =  gcr( SNR, GAMMA) 

%  This  function  computes  the  collapsing  ratio  caused  by 
%  insufficient  video  bandwidth  in  a  Gaussian  square  law 
%  system.  GAMMA  can  be  a  vector  but  SNR  must  be  single 
valued. 

%  Note  that  GRHO  is  also  the  equivalent  number  of  independent 
%  samples  averaged  by  the  video  filter. 

%  Reference:  Thesis  Eqn  (26). 

% - 

B12  =  1  +  2*GAMMA  .^2 
Bll  =  1  +  GAMMA  .^2 

B  =  2*SNR*sqrt(B12  ./BH)  ”  sqrt(B12) 

O  =  -2*SNR*sqrt(B12) 

GRHO  =  0.5*(-b  +  sqrt(B  .^2  -  4*0) 


7.  HMT.N  FUNOTION 
function  H  =  hmt(Z,m) 

%  This  function  computes  the  mth  order  Heirmite  polynomial  for 
%  values  of  m  less  than  or  equal  to  12. 

% - 

if  m  ==  0 
H  =  1 

elseif  m  ==  1 
H  =  Z 
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2 


elseif  m  == 


H  =  Z  .^2 

- 

1 

elseif  m  == 

3 

H  =  Z  .^3 

- 

3*Z 

elseif  m  == 

4 

H  =  Z  .^4 

- 

6*Z  . 

■  3 

elseif  iti  == 

5 

H  =  Z  .^5 

- 

10*Z 

.  "3 

+ 

15*Z 

elseif  m  == 

6 

H  =  Z  .^6 

- 

15*Z 

.  *4 

+ 

45*Z  .^2  -  15 

elseif  m  == 

7 

H  =  Z  .^7 

- 

21*Z 

.  "5 

+ 

105*Z  .^3  -  105*Z 

elseif  in  == 

8 

H  =  Z  .^8 

- 

28*Z 

.^6 

+ 

210*Z  .^4  -  420*Z  .^2  + 

elseif  m  == 

9 

H  =  Z  .^9 

- 

36*Z 

.  ^7 

+ 

378*Z  ."5  -  1260*Z  .^3 

elseif  m  == 

10 

H  =  Z  .^10 

- 

45*Z 

.  "8 

+ 

630*Z  .-'6  -  3150*Z  .^4  ^ 

-  945 

elseif  m  == 

11 

H  =  Z  .'^11 

- 

55*Z 

.  ^9 

+ 

990*Z  .^7  -  6930*Z  .''5  4 

-  10395*Z 
elseif  m  ==  12 

H  =  Z  .^12  -  66*Z  .^10  +  1485*Z  .^8  -  13860*Z 
+  51975*Z  .^4  -  62370*Z  .''2  +  10395 

end 


105 

+  945*Z 
-  4725*Z  ."2 

■  17325*Z  .''3 

.^6 
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